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At large virtuality Q 2 , the coupling to the vector meson production channels provides us with 
a natural explanation of the surprisingly large cross section of the neutral pion electroproduction 
recently measured at Jefferson Laboratory, without destroying the good agreement between the 
Regge pole model and the data at the real photon point. Elastic rescattering of the n° provides 
us with a way to explain why the node, that appears at t ~ —0.5 GeV 2 at the real photon point, 
disappears as soon as Q 2 differs from zero. 
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The electroproduction of neutral pion raises two is- 
sues. The first one is a long standing issue (mid 70's): 
the node that appears at t ~ -0.5 GeV 2 at the real pho- 
ton point l| disappears as soon as the virtuality of the 
photon Q 2 differs from zero 0. Several attempts, rang- 
ing from Regge cuts [l-Q to direct coupling to quarks 
and rearrangement Q , were proposed to qualitatively ex- 
plain this result. However, no quantitative explanation 
has been proposed yet, except at the expense of a strong 
variation with Q 2 of the elatic cut 0] • The second issue is 
more recent and was raised, about three years ago, during 
the preliminary analysis of the data obtained at Jefferson 
Laboratory (JLab) in HallA and in HallB. When extrap- 
olated at large Q 2 , any Regge pole model that leads to a 
fair understanding of the low Q 2 data grossly underesti- 
mates the JLab data. The HallA data [8j have just been 
released; the HallB data @ are still in the final stage of 
analysis and will be released soon. 

In this Letter, I propose a solution of these two prob- 
lems which combines elastic ir° rescattering and the cou- 
pling to inelastic up, p + n, p + A and p~ A ++ channels. A 
recent measurement at JLab [Io| has shown that the cross 
section of the p("f* , p + )n channel, which is very small at 
the real photon point, becomes comparable to, and even 
larger than, the cross section of the p(j*,p°)p channel, 
at large virtuality Q 2 . 

Figure [1] summarizes these findings. The DESY data 
have been recorded in a geometry where the ir° is emit- 
ted perpendicular to the electron scattering plane (az- 
imuthal angle (j> — 90°) that emphasizes the transverse 
component of the cross section. The JLab data have been 
integrated over the tt° azimuthal angle and correspond to 
the combination of the transverse and longitudinal cross 
sections dax/dt + tda^jdt, where e is the polarization 
of the virtual photon. The combination of the Regge 
pole and the tt° elastic scattering amplitude (Pomeron 
cut) allows to reproduce the cross section at moderate 
Q 2 (~ 0.5 -i- 1 GeV 2 ), but underestimates the JLab data 
at higher Q 2 . The coupling to the charged p production 
channels brings the theoretical cross section close to the 
JLab data at Q 2 = 2.3 GeV 2 . It is remarkable, and was 
rather unexpected, that the cross section, at — t ~ 0.3, 



does not change between Q 2 — 0.85 and 2.3 GeV 2 . 
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FIG. 1: (Color on line) The cross section of the p(7*,7r°)p 
reaction recorded at DESY top and bottom left, and 
JLab [H, bottom right. The basic Regge Pole model with the 
Pomeron cut corresponds to the red (without 6i pole) and the 
black (with 6i pole) full line curves. The dashed line curves 
(green) take into account the contribution of the n charge ex- 
change scattering cuts, the inelastic up and p + n cuts. The 
(blue) full line curves take also into account the contribution 
of the inelastic p*A cuts. 

The basic Regge pole amplitudes are fully described in 
ref. [IH (GLV). Instead of their expression in terms of 7 
matrices, I use the expression of the u and p (that are 
given in the appendix of ref. 0), as well as the expres- 
sion of the b\, t-channel exchange amplitudes in terms 
of a matrices. In the GLV scheme, the node in the real 
photon cross section near t = -0.5 GeV 2 is naturally gen- 
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erated by the use of a non-degenerated amplitude of the 
u Regge pole. The contribution of the p Regge pole with 
a degenerated amplitude fills in the dip and brings the 
GLV model close to the data. The exchange of the b\ me- 
son does not contribute significantly to the unpolarized 
cross section, but is needed to account for the measured 
photon asymmetry £. Multiplying each amplitude by an 
electromagnetic form factor Fm(Q 2 ) keeps the shape of 
the real photon cross section in the virtual photon sector, 
in contradiction with experiment. 

The other way to generate a node in a Regge ampli- 
tude is to use a degenerated amplitude and supplement it 
by the elastic absorptive cut. Assuming that the elastic 
scattering amplitude is driven by the Pomeron exchange, 
it is possible to express the cut amplitude as an effective 
Regge pole [14J. In this scheme, the u Regge exchange 



amplitude takes the form: 
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where t={k u 



7 ) 2 is the four momentum transfer and 
m w is the mass of the oj meson. The Feynman amplitude 
Tf has the same spin-momentum structure as in the GLV 
scheme. Only the g U NN coupling constant will be re- 
fitted to experiment. The u> Regge trajectory is the same 
as in the GLV scheme: 



,(0) 



t = 0.44 + 0.9t 



(2) 



while the intercept and the slope of the effective trajec- 
tory of the cut take the form: 



a c (0) = a u (0) +a P (0) - 1 =0.44 
< = (a' u xa' P )/(al + a' P )= 0.2 



(3) 



where the intercept and the slope of the Pomeron Regge 
trajectory are respectively ap(0) — 1 and a' p = 0.25 

The purely destructive interference between the pole 
and the cut amplitudes is a direct consequence of the 
structure of the rescattering loop amplitude and the al- 
most purely absorptive nature of the elastic wN scatter- 
ing amplitude (see eqs. (9) to (11) and the corresponding 
discussion in ref. [151]). It is possible to reproduce the 
GLV non-degenerated u amplitude (squared) with the 
following choice (dotted line curve in Fig. [2]): 



G c (t) 
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It is worth pointing out that the coNN coupling constant 
is about half of the GLV one {g^NN/^ = 17.9), in better 
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FIG. 2: (Color on line) The cross section of the p(j, ir°)p reac- 
tion [lj (top) and the photon asymmetry (bottom). The 
dashed line curve corresponds to the uj Regge Pole amplitude 
with a degenerated trajectory. The dotted line curve includes 
the Pomeron cut. The red dot-dashed line curves include the 
p exchange, while the black full lines curves include also the 
b\ exchange. The blue short-dashed line curves include all the 
inelastic cuts. 



agreement with the range of values that are determined in 
the analysis of low energy NN scattering data. It is also 
almost the same as the value (ff 2 ; jvjv/^ :7r = 6.44) needed in 
the analysis of the p(7, rf)p reaction (l7| . In this channel, 
there is no node and the r]N scattering cross section is 
much lower than the irN one: The use of a degenerated 
u> Regge amplitude alone, with no elastic cut, is more 
justified. The exchange of the p degenerated trajectory, 
with the same coupling constants as in GLV, brings the 
model close to the unpolarized cross section. 

In the first version of this paper 18( an axial- vector 
biNN coupling was used, as in the GLV scheme. It turns 
out that this coupling is forbiden: it is allowed at the 
aiNN vertex, but the a-i meson does not couple to the 
neutral pion. Therefore, I use the following b\ exchange 
current which is based on an axial-tensor biNN coupling: 



■?bi = ^^9b 1 ^u(pf)^ 5 u(p l )(2p> 1 f kp l -2pi*kp f ) 



-iira^ (t) 



so 
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where the Regge trajectory is the same as the pion degen- 
erated trajectory and the electromagnetic coupling con- 
stant <7b l7 T7 is the same as in GLV. The choice of the 
strong coupling constants of the b\ meson g^/Aw =21.5 
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and Kb 1 = 2 leads to a fair agreement with the beam 
asymmetry at the real photon point. 

This scheme offers us with a way of shifting the min- 
imum of the cross section when the virtuality Q 2 of 
the photon increases, by using slightly different cut off 
masses in the electromagnetic form factors of the poles 
and the Pomeron cut, F(Q 2 ) = 1/(1 + Q 2 /A 2 ). The fol- 
lowing choice leads to a good accounting of the DESY 
data at Q 2 = 0.85 GeV 2 : A 2 = 0.325 GeV 2 , A 2 = 
0.400 GeV 2 , A 2 bi = 1. GeV 2 , A 2 = 0.300 GeV 2 and 
S C (Q 2 ) = -0.46g 2 /0.85. 

The agreement is good too at Q 2 — 0.55 GeV 2 , but 
it is not possible to get rid of the second maximum in t 
when Q 2 — 0.22 GeV 2 . Also, the extrapolation of this 
scheme at Q 2 = 2.3 GeV 2 misses the recent JLab data. 
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FIG. 3: The Charge Exchange pion rescattering graphs. 



The first inelastic cut that may play a role is the 
Charge Exchange (CEX) pion rescattering cut (Fig. [3]). 
The charged pion electroproduction [ujj . around Q 2 = 
2.3 GeV 2 , is larger (~ 3 pb/GeV 2 ) than the tt° one (~ 0.4 
pb /GeV 2 ) at low t. Neglecting its principal part, the cor- 
responding rescattering matrix element reduces to [lij ]: 
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where p' cm = y/ (s — (m ff — m) 2 )(s — (m^ + to) 2 )/4s is 
the on-shell momentum of the intermediate neutron, for 
the cm. energy \fs. The two fold integral runs over 
the solid angle Vt of the intermediate neutron, and is 
performed numerically. The four momentum transfer 
between the incoming photon and the intermediate 7r 
is t 1 = (kry — -Ptt) 2 , while the four momentum trans- 
fer between the intermediate and the outgoing pions is 
t„ = — -Pn-) 2 - The summation over all the spin indices 
of the intermediate particles is meant. 

For the p(-f* ,tt + )ti amplitude, I use the VGL 
model 20] which reproduces fairly well the experimen- 
tal datapj] around Q 2 =2.3 GeV 2 and y/s = 2.2 GeV, 
at least the Longitudinal part. The expression of the 
CEX amplitude is: 



Tcex 



3g p (l + k v ) 



(7) 



where g 2 T7r /47r = 5.71, where g 2 /47r = 0.92, where Ky = 
6 and where F\ is the nucleon form factor as defined 
in [15j. Since the experimental t distribution exhibits a 
node, I use the non degenerated Regge propagator V P R 
with the saturating trajectory of the p |12[ . As shown in 
Fig. |H this gives a good account of the 7rN CEX scatter- 
ing data 21( in the same energy range. 
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FIG. 4: (Color on line) The ttN CEX at s/s = 2.4 GeV [21[. 
The (black) line, marked "lin", corresponds to the use of a 
linear Regge trajectory, while the (red) line, marked "sat", 
corresponds to the use of a saturating Regge trajectory. 



Since the ttNA and pNA coupling constants are com- 
parable to and even larger than the irNN and the 
pNN ones (see for instance ref. 22[), the 7r + A° and 
7r~A ++ intermediate states play also a role (Fig. [3]). The 



rescattering matrix element is a straightforward exten- 
sion of cqs. (H|) and ([7|), using the relevant coupling con- 
stants as well as the relevant isospin coefficients, and re- 
placing the a matrices by the N — ¥ A spin transition 
matrices S. 

Fig. [5] shows that the coupling to these CEX channels 
is not enough to account for the large experimental cross 
section at Q 2 = 2.3 GeV. 

The next cuts are the vector meson cuts shown in 
Fig. [6] The generic amplitude is: 
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where p c . m . = y/(s - (m v - M) 2 )(s - (my + M) 2 )/4s 
is the on-shell momentum of the intermediate baryon (of 
mass M), for the cm. energy s/s. The two fold integral 
runs over the solid angle f2 of the intermediate baryon. 
The four momentum transfer between the incoming pho- 
ton and the vector meson is i 7 = (fey — Py) 2 , while the 
four momentum transfer between the vector meson and 
the outgoing pion is t n = (k^ — Py) 2 ■ The summation 
over all the spin indices of the intermediate particles is 
meant. 
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FIG. 5: (Color on line) The contribution of the ttN (dashed 
red line), the 7rA (full black line) CEX cuts and of the cop cut 
(dash-dotted black line) to the cross section at Q 2 = 2.3 GeV 2 
and yfs — 2.269 GeV. The dotted blue line is the contribution 
of the Regge poles. The contribution of the ttN and the tiN 
+ 7rA cuts alone is shown in the bottom of the figure. The 
dotted lines correspond to the longitudinal component only. 



where / 2 /47r = 18.4, and where the amplitude is eval- 
uated with the actual kinematics of the p(ui,ir )p reac- 
tion. This model leads to a very good agreement with 
the available data (Fig. [7]). Again, the contribution of 
the uip cut is not enough to reproduce the experimental 
data at Q 2 = 2.3 GeV 2 (Fig. [5]). 




FIG. 7: The cross section of the p(ir , uj)n reaction at ^/s — 
2.785 GeV HJ. 
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FIG. 6: The Vector meson cuts. 
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Since the p° cannot decay into two tt s, only the tup, 
the p + n and the p ± A cuts have to be taken into account. 
In the ujp cut, the amplitude of the p(j* , to)p reaction is 
based on the exchange of the Regge trajectories of the 
Pomeron, the 7r, the ji in the t-channel and of the pro- 
ton in the u-channel. The model is described in ref. [23I ] 
and reproduces well the experimental data 24| in the 
JLab energy and momentum range. The amplitude of 
p(ui,ir°)p has the same structure [l3j as the p exchange 
part of the Regge amplitude of the reaction ^(7, ir°)p, to 
which it is related under the Vector Meson Dominance 
assumption: 



fu, 
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FIG. 8: (Color on line) The comparison between the cross 
sections of the p(7, p°)p (black) andp(7,p + )n (red) reactions 
at Q 2 = (top) and Q 2 = 2.3 GeV 2 (bottom). The open red 
circles and the dash-dotted red curve correspond to the photo 
production of p~ , while the filled red circles and the full red 
curves correspond to the electroproduction of p + . 



The contribution of the p + n cut is far more impor- 
tant (Fig.[T|). The first reason is that the cross section of 
the N(p. tt)N reaction is larger than the cross section of 



the N(uj,ir)N reaction at low t (compare Fig. 3 of 15] 
and Fig. [7]): The former reaction is driven by tt exchange 
while the latter is driven by p exchange. The second rea- 
son is that the cross section of the p(-f, p + )n, which is very 
small at the real photon point, becomes large at large Q 2 . 
This is shown in Fig. [8] which compares the cross sections 
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FIG. 9: The graphs of the reaction p(j, p + )n.Top: meson 
exchange graphs. Bottom: Contact and nucleon exchange 
gauge graphs. 



of the p(7, p°)p and p(j, p + )n reactions at the real pho- 
ton point and at Q 2 = 2.3 GeV 2 . The reference to the 
experimental cross sections of the p°, as well as the de- 
scription of the model, can be found in [23], [2(| . At the 
real photon point, the cross section of the 71(7, p~)p reac- 
tion comes from |27|. At Q 2 = 2.3 GeV 2 , the two data of 
the p(7, p + )n reaction are preliminary and given for illus- 



tration only. They come from the preliminary analysis of 
a JLab experiment (lp| . The model of the 73(7, p + )n re- 
action is based on the exchange of the 7r + and p + Regge 
trajectories. The 7r exchange amplitude is the same as 
the amplitude [23, 26| that dominates the p(7*, uj)p reac- 
tion (with trivial changes of the coupling constants). It 
is not important and the dominant contribution comes 
from the p exchange amplitude. Its vector part takes the 
form: 
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where e and f p are respectively the polarization vectors 
of the incoming photon and the outgoing p, and where 
the coupling constants are g 2 /47r = 0.4 and ny = 6. One 
recognizes easily the p pole term as well as the nucleon s- 
channel pole term and the contact term (see Fig. [9]) that 
are necessary to make the amplitude gauge invariant. 

I use the GLV degenerated Regge propagator V P R , with 
a linear trajectory, and the t dependent electromagnetic 



form factor F p (Q 2 ,t) of 23[ with a cut-off mass of A 2 



0.9 GeV 2 . Under those assumptions, the model not only 
predicts the right transverse and longitudinal integrated 
cross sections (Fig. [5]) from the real photon point to 
large Q 2 , but also their t distribution in the whole range 
(1.5 < Q 2 < 4 GeV 2 , 2 < Vs < 2.8 GeVjthat has been 
covered in the JLab HallB experiment [10(. The model 



will be compared to the final JLab data in the experi- 
mental paper. 

As shown in Fig. [TJ the contribution of the p + n cut 
accounts for about half of the cross section at Q 2 = 
2.3 GeV 2 , modifies little the cross section at Q 2 = 0.55 
and 0.85 GeV 2 , but interferes with the Regge pole con- 
tribution at Q 2 — 0.25 GeV 2 . At this point of the discus- 
sion it is worth to emphasize that no freedom is left in the 
amplitudes of the cuts, in that sense that the unitarity in- 
tegral relates on-mass shell elementary amplitudes of the 
production and absorption of the intermediate mesons. 
As soon as those amplitudes reproduce the elementary re- 
action cross sections, the amplitude of each cut is frozen 
and their relative contribution is driven by the actual 
relative size of the cross section of the elementary chan- 
nels. So, this singular part of the rescattering integrals is 
on solid grounds and consistently relates several channels 
that have been studied recently at JLab or elsewhere: it 
simply cannot be overlooked. 

The contribution of the principal part of the rescat- 
tering integrals is not so well constrained as it requires 
the knowledge of either the off-mass shell amplitudes, if 
is calculated by brute force, or their asymptotic behav- 
ior, if it is determined by a dispersion relation. This is 
an open issue, but I simply note that the measured ratio 
between the teal part and the imaginary part of the am- 
plitude does not exceed 10% in Compton scattering or 
meson-Nucleon elastic scattering at forward angles. So 
there is good reason to expect that the principal part of 
the rescattering integrals contributes little. 

The third reason why the coupling to the charged p 
production channels is strong is that the contribution of 
the A's intermediate states (Fig. [6]) is as important as 
the contribution of the neutron intermediate state. Since 
there are no measured cross sections of the 7V(7*,p ± )A 
reaction in the virtual photon sector, I note that, at low- 
est order, the structure of the lower part of the loop di- 
agrams is very similar: 

(7 ' kpcr ' — kp ■ A>7r H \~ i(7 ' kp x kjr 

' kpS ' /C-tj- = 77^p ' ^7T ^~n^ ' ^ ^ 7r (H) 

Where S is the spin operator of the N — > A transition. 
Assuming that the scalar and vector parts contribute 
equally: 
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where p A and p are the momenta of the A and the neu- 
tron in the pA and p + n loop respectively, and where the 
ratio of the coupling constant is G p G v / 'g p g v — 4.49 ac- 
cording to [22j . The last bracket contains the ratio of the 
isospin coefficients (1/2 in the A ++ channel and 1/6 in 
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the A channel). For the 7rA cuts, such an estimate is 
consistent with the numerical evaluation of the integral 
that is shown in the bottom part of Fig. [5] I also note 
that, at the real photon point, the same kind of estimate 
predicts a cross section of the p(7,p + )A° reaction com- 
parable to the cross section of the ivOy, p + )n reaction, in 
accord with the only experiment [281 ] available so far. 

Under those assumptions, the model becomes close to 
the JLab HallA data at Q 2 = 2.3 GeV 2 (Fig. [I) without 
destoying the good agreement at the real photon point 
(Fig. [2]). It slightly overestimates the DESY data at Q 2 = 
0.55 and 0.85 GeV 2 . This can be fixed by a fine tuning 
of the cut off masses in the electromagnetic form factors 
of the Regge poles, which I had chosen to reproduce the 
DESY data without the cut contribution. This will not 
affect the agreement at the lowest and the highest Q 2 . 
However the model already gives a good account of the 
preliminary JLab HallB data |9| in the range 1.5 < Q 2 < 
4 Gev 2 , 2.1 < y/s < 2.8 GeV, and I postpone such an 
adjustment until the final data are released. 
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FIG. 10: (Color on line) The response functions of the reac- 
tion p(7*,7T°)p in the two JLab HallA kinematics settings. 
Dashed lines: pole contributions and Pomeron cut alone. 
Dash-dotted lines: without pA cuts. Full lines: pA cuts in- 
cluded. 



2tt- 



da 

dtdcj) 



dt 



+ecos(2(j))aTT + \/2e(e + 1)cos{4>)(Jtl 
+hy/2e(l - e)sin{(j))a' TI }j (13) 

and I have renormalized the experimental tranverse- 
longitudinal cross sections [§] , <jtl ad a' TL , accordingdly. 

Besides the unpolarized cross section (<tt + eox), the 
model reproduces also the Transverse- Transverse inter- 
ference response function ott- The understanding of 
the slow variation of the cross section is a non triv- 
ial result. At the lowest Q 2 = 1.94 GeV 2 , the mass 
of the intermediate state is y/s — 2.052 GeV, slightly 
above the pA threshold (2.010 GeV). At fixed X = 
0.368, y/s increases when Q 2 increases, and the onset 
of the pA cuts compensates the decrease with Q 2 and 
y/s of the Regge amplitudes that drive the loop ampli- 
tudes. The contribution of the cuts is small in cttt, 
but large in the Transverse-Longitudinal response func- 
tions otl and cr' TL . The model predicts the same sign 
and magnitude as the JLab HallA fifth response function 
a' TL (polarized electrons), and therefore reproduces the 
beam asymmetry Alu (Fig. E]) that has been recorded 
in HallB (25|. But it predicts a different sign for the 
Transverse-Longitudinal response function cttl- Since it 
predicts the same sign and magnitude as the preliminary 
HallB utl i, I leave open the discussion until the final 
data are released. 
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FIG. 11: The beam asymmetry 



recently measured at 



JLab, for Q 2 = 2.3 GeV , X = 0.35 (y^ = 2.27 GeV) and , 
90°. The dotted curve takes only into account the uj, p and 
bi Regge pole contributions. The dash-dotted curve includes 
also the CEX cuts contribution. The full curve includes all 
the cuts contributions. 



Fig. Uni shows the comparison of the model with the 
response functions which have been determined at JLab 



Finally, the unitarity cuts are also the key to the 
successful interpretation of the target asymmetry Ajjl 
(Fig. [T2j!. The CEX cuts contribute, but falls short. The 
p + n and pA cuts lead to the right sign and the right 
magnitude of the signal. 
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FIG. 12: The target asymmetry 30j recently measured at 
JLab, for Q 2 = 1.82 GeV 2 , s/s = 2.358 GeV and t = - 
0.31 GeV 2 . Dashed line: lj and p poles. Dotted line: 61 
pole included. Dash-dotted line: tt CEX cuts included. Full 
line: Vector meson cuts included. 



In conclusion, the coupling to the vector meson pro- 
duction channels provides us with a natural explanation 
of the large cross section of the ir° electroproduction 
cross section at Q 2 ~ 3 GeV 2 . Only a few intermediate 
hadronic states contribute to the unitarity cuts, which 
are on solid grounds when the elementary production 
and absorption cross sections are large and are known. 
So far, we are still in the hadronic regime and we have 
not reached the domain of factorization between a hard 
perturbative scattering and a soft non perturbative nu- 
cleon structure function. 
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